February 1, 2008 5:1 WSPC/INSTRUCTION FILE Proc 



Modern Physics Letters A 

© World Scientific Publishing Company 



in 
o 

O ' MONTE CARLO HAMILTONIAN OF LATTICE GAUGE THEORY 

F. PARADIS 

, Department de Physique, Universite Laval, 

Quebec, Quebec, G1K 7P4, Canada 
email: fparadis(3phy.ulaval.ca 

H. KROGER* 

Department de Physique, Universite Laval, 
■ Quebec, Quebec, G1K 7P4, Canada 

email: hkroger@phy.ulaval.ca 

o 

oo 
O 

\f} , Department of Physics, Zhongshan University, 

i Guangzhou 510275, People's Republic of China 

email: stslxq@zsu.edu.cn 

K.J.M. MORIARTY 

Institute for Computational Studies, 
Department of Mathematics, Statistics and Computational Science 
' Dalhousie University, Halifax, Nova Scotia B3H 3J5, Canada 

email: kevin.moriarty@cs.dal.ca 



- r— I 

X 

5-S 



X.Q. LUO 



Received (Day Month Year) 
■ Revised (Day Month Year) 



We discuss how the concept of the Monte Carlo Hamiltonian can be applied to lattice 
gauge theories. 

Keywords: Lattice gauge theory; Hamiltonian formulation. 
PACS Nos.: 11.15.-q, 03.70.+k 



1. Bargman states 

In standard lattice field theory it is customary to compute vacuum expectation 
values of operators, like 

< Sl,t = +oo\A(x,t)B(x',t')\n,t = -oo > , (1) 

which denotes the amplitude that the system is in its ground state |f2 > at t = — oo, 
is acted upon by some local field operators at some intermediate time and returns to 
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the ground state at t = +oo. In lattice calculations, using a lattice of finite extension 
in space and time, one can extract from such correlators information about particle 
masses in the ground state and in the lowest lying excited states. In order to extract 
information about particle properties in higher excited states one may consider the 
following amplitudes 

<* fi ,t = T\* in ,t = 0> , (2) 

where |\& > denotes a Bargman state. A Bargman state on a space-time lattice at 
a given time slice is defined by assigning a particular value of a local field variable 
to each lattice node. In the case of the real scalar bosonic field $>(x, t), for instance, 
in a given time-slice t each lattice node Xk is associated with some real number 
(taking values between +oo to — oo). This may be visualized as a regular array in 
a plane with perpendicular spikes (of different length) poking out of the plane on 
either side. The ground state can be considered as a linear superposition of many 
Bargman states. On the other hand, a single Bargman state can be considered as 
a linear superposition of the ground state plus of many excited states. Thus, by 
computing matrix elements like in Eq.JSJ, one expects to get information on the 
sector of excited states of the system. The reason why we have chosen Bargman 
states in particular, is that those are localized states, which is very convenient in 
local quantum field theory. 



2. How to compute transition amplitudes by Monte Carlo? 

The idea is to express the quantum transition amplitude in terms of a path integral 
and to evaluate such path integral numerically, using Monte Carlo with importance 
sampling. However, such Monte Carlo method allows only to compute the ratio of 
transition amplitudes, like 

R = <*fi,t = T\6{t = T int )\V in ,t = > 

but not a single transition amplitude. A solution of this problem has been proposed 
for quantum mechanical system^, and also for the scalar model in field theory^. 
The idea is to split the action occuring in the path integral into a free, i.e. kinetic 
term plus an interacting, i.e., potential term. As an example let us consider a 1-d 
quantum system of a massive particle in some potential 

S[x}= [ dt\mx 2 - V(x) = [ dt\mx 2 - ( dtV(x) = S kin [x] + S pot [x] .(4) 
Jo 2 J 2 J 

Then the transition matrix element in imaginary time can be expressed (we use 
H = 1 throughout) 



< z\exp[-HT]\y > = < z\exp[-H km T]\y > 

J[dx] exp[-S pot [x}] exp[ 



|2,T 



J[dx] eM-S km [x]]\ Z y T 



(5) 
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Here exp[— S pot ] is treated as an observable and exp[— S ] represents a weight fac- 
tor. The ratio of path integrals can be computed by standard Monte Carlo methods 
with importance sampling. The matrix element < z\exp[— H km T]\y >, describing 
the free propagator (in imaginary time) is known analytically, 



< z\exp[-H« m T]\y >= x j — exp - — (z - yf . (6) 



m 



2T 



The calculation of this transition amplitude proceeds most simply by using momen- 
tum space, which diagonalises H km given by 

p2 

2m 

The connection between position space and momentum space is given by 

< x\p > = exp [ipx] . (8) 

\/27r 

Note that this scalar product is solely determined by the algebraic commutator 
relation of the operators X and P, plus the normalisation of the position and mo- 
mentum states. 



3. Splitting the action in lattice gauge theory 

Let us consider QED on the lattice without fermions. The group of gauge symmetry 
is U(l). We keep in mind how to construct a lattice Hamiltonian (see e.g, Creuta^l 
Rothe^l) via the transfer matrix, by splitting the lattice action into terms of time- like 
and space-like plaquettes, 

SP} = \- £ [l-i*e(f/ F )] + 4- E [l-Jfe(tfp)] 
5 fl o p 9 a p 

^time-like ^space-like 

= S km [U} + S pot [U] . (9) 



The transition amplitude can be split up, in analogy to the example from quantum 
mechanics, Eq. |JSJ, 



< U fi \exp[-HT]\U in > = < U ft \eM-H km T]\U m > 

jkin\ T m\ U f^ T 

(10) 



J[dU] exp[-S pot [U]] exp[-S kin [U]]\ v (1 



J[dU] exp[-S kin [U}}} 



Ufi,T 
U in ,Q 



Like in the quantum mechanical example, the ratio of path integrals can be com- 
puted using Monte Carlo with importance sampling. The problem is now the com- 
putation of the transition amplitude evolving according to the kinetic term. 
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4. Transition amplitude of evolution due to kinetic term 

In order to calculate the amplitude 

< U ft \exp[-H km T]\U m > 



(11) 



we make use of the close correspondence between quantum mechanics and lattice 
gauge theory. There is a lattice Hamiltonian, which is the generator of infinitesimal 
translations in time, the so-called Kogut-Susskind Hamiltonian. In the case of U(l) 
gauge theory, it reads ^ 



(12) 



Here Uj denotes the operator of the electric field between spatial neighbour lattice 
sites i and j . This operator on the lattice has discrete eigenvalues, 



hjWj >— MjWj > j — 0j il; = t2, . . 
We can express the transition amplitude, Eq. iJTTJ, by 

<U fi \ex.p[-H kin T]\U in > 



(13) 



= J2< U fi \\ >< A| exp[~H km T}\\ >< X\U in > 



< Ufi\X >< A|exp 



A, A' 



9 2 T 



EH 



i 2 - 

2a ^ 13 

<ij> 



\X' >< X'\U in > 



n e 



< C^|Ay > exp 



=0.±1,±2.. 



9 2 T 2 



<\ l3 \U™>> ■ (14) 



Due to the locality of the Bargman-link state and the structure of the kinetic part 
of the Kogut-Susskind Hamilton operator, this amplitude factorizes in separate 
amplitudes for each link. In order to evalute such individual link amplitudes one 
needs to compute the scalar product which defines the transition from the link basis 
to the electric field basis, 



< UijlXij > 



(15) 



4.1. Scalar product < C7|A > 

Let us compute the scalar product given by Eq. (|15l) . For simplicity of notation, we 
drop the subscript ij in this section. Like in quantum mechanics, this scalar product 
is determined solely by the algebraic commutator relation between the operators U 
and I, plus the normalisation of the link and electric field states. Let us recall the 
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properties of the link basis 



3 



u\u >= u\u > 

< U'\U >= 5{U' - U) 

1 = J dU\U><U\, (16) 

where J dU denotes a group integral and S(U' — U) denotes a 5 function (distribu- 
tion) over group elements. Likewise, the electric field basis obeys 

l\X >= A|A > 
< A'|A >= 6 x , lX 

1= |A><A|. (17) 

A=0,±l,±2,... 

The algebraic commutator relation of 1 and U is given b}GE] 

[l,U]=-U. (18) 
First, we want to compute the matrix element < U'\l\U >. Using Eq. (|18|) . we have 

< u'\[i,u]\u >=< u'\iu\u > - < u'\ui\u >= [u-u'\ < u'\i\u > , (19) 

and on the other hand 

< U'\[l,U]\U >=< U'\ - U\U >= -U < U'\U >= -U8{U' - U) . (20) 
Thus we find 

mv - m b 

(21) 



< u'\i\u >= -u s{ v' g = Hu^tjW u) 



U'-U x '~ dU 
where we used the following result from distribution theory 



8{x' - x) 



(22) 



x' — x dx' 

Second, we want to compute the matrix element < A|£7 >. Using Eq. (|18f) . we obtain 



and also 



< X\[l,U]\U >=< A| - U\U >=-U< X\U > , 

< x\[i,u]\u >= xu < x\u > — < x\ui\u > 
= xu < x\u > - I dV < x\u\v >< v\i\u > 



(23) 



= XU < X\U > + J dV UV < X\V > -^&{V - U) 
= xu<x\u>-^ v [UV<X\U>]\ V=U 



xu < x\u > 



u <x\u> +u 2 — < x\u > 

oU 



(24) 
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Combining Eqs. (|23I24[1 yields the differential equation 

U^j < X\U >= X< X\U > . (25) 
The solution has the form 

< X\U>=B X (U) X , (26) 

where B\ is some constant. This constant is determined from the normalisation 
condition of the electric field states, Eq. (|17|l 

5\i \ =< A'|A >= J dU < X'\U >< U\X > 

= By (B x )* J dU (U) x '- X = By (Bx)*-^ J da exp[ia(X' - A)] 

= By(B X )*5y,X , (27) 
Here we have used the orthogonality of periodic functions 

^ e ine e -in'C = 2ir8 n , n , . (28) 



We have parametrized a link U = exp(ia), with the group integral measure (Haar 
measure) dU = da/(2n). From Eq. (|27|l we conclude 

B x = l. (29) 

From Eqs. (|26I29|I we conclude 

< X\U >= (U) x . (30) 
4.2. Transition amplitude for a single link 

Now we obtain the transition amplitude for a single link, denoting Ufi = exp(za/i) 
and U ln = exp(«a m ) 

<U fi \ exp[-H kin T]\U in > 

<U fi \X>e^[~^X 2 ]<X\U m > 

A=0,±l,±2,... 

g 2 T 

= ^2 ex P[ — ^-A 2 ] exp[iA(a in - a/i)] 

A=0,±l : ±2,... 

= ^2 exp[ — ^-^ 2 ] cos[n(ai„ - a fi )} . (31) 

n=0,±l,±2,... 



Let us do a consistency check of this result. Let us consider the time evolution 
during a time interval T — ao (lattice spacing in time) and then take the limit 
a -> 0. 
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We define 



Aa = 



A 2 = 9 2 ao 



2a 

= An, n = 0,±1,±2, Ax = A . (32) 



Then the transition amplitude, Eq. l|31fl becomes 
< Ufi\ exp[-H km ao}\U in > = — ^ Aa; exp[-x£] cos[x„Aa/Aa 



Ax n=u. i. 2 



1 / >+ °° 

— / da; exp[— x 2 ] cos(xAa/A) 

A J-oo 



/tt (Aa) 



A^ AA 2 



2ira a (Aa) 



2 ex P 2 o J 

~oo--o 2n8(afi - a in ) = S(U fi - U in ) , (33) 

which reproduces the correct behavior in the limit when ao goes to zero. Moreover, 
in the second to last line one may replace (a/j — ai„) 2 /2 by 1 — cos[a/i — a^]. 
This is justified in distribution theory when ao goes to zero. Then one obtains the 
expression of the path integral of link variables for a single time step. 

It is interesting to see an application of this approach to non-abelian lattice 
gauge t heory and dynamical quarks. Hamiltonian lattice QCD at finite chemical 
potentia l I I I might be a more natural solution to the infamous complex action 
problem. 
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